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1. INTRODUCTION 



Consider the fourth-order ordinary differential equation 



a; e R, 



(1.1) 



where the prime denotes the derivative with respect to the independent variable and 
the potentials u and v are real valued in such a way that w, tt', v are intcgrable and they 
decay exponentially or faster [8,9] as a; — > ±oo. Let us write (1.1) as C(p = /c^t/?, where C 
is the linear operator defined as 



We will consider (1.1) in the sector f2 in the complex A;-plane where argA; e [0,7r/4] with 
the understanding that we exclude the point /c = 0. A complete analysis at A; = for (1.1) 
does not exist and this deserves a careful future study. Some partial analysis at /c = is 
available in [8] . We will omit the analysis at A; = in this paper. 

The differential equation appearing in (1.1) is the canonical equivalent [4,6] of the 
Euler-BernouUi equation, and the former is obtained from the latter via a transforma- 
tion of both independent and dependent variables [4] . Such fourth-order equations arise 
in analyzing vibrations of beams, whereas second-order equations are used in describing 
vibrations of strings. 

We refer the reader to [8,9] and [7,10] for two complementary studies of (1.1). Under 
appropriate restrictions on the potentials u and f , Iwasaki [8,9] studied (1.1) by analyzing 
it in the sector VL under the assumption that there are no bound states and no spectral 
or nonspectral singularities. He obtained [8] various properties of solutions to (1.1) and 
formulated [9] the inverse problem of recovery of u and v from some appropriate scatter- 
ing data. Iwasaki posed [9] the inverse scattering problem for (1.1) as a boundary- value 
problem where the scattering data consisted of a reflection coefficient and a connection 
coefficient specified on the boundary of and he provided the proof of uniqueness for 



C := 



D"^ + DuD + V, 



dx 



(1.2) 
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the solution to that inverse problem. A special case of (1.1) was studied in [7,10] under 
the very restrictive assumption that the reflection and connection coefficients are all zero. 
The corresponding inverse problem was formulated [10] as a Riemann-Hilbert problem 
on the whole complex /c-plane, where the set of scattering data is specified on the rays 
arg A; = {I — l)7r/4 for I = 1, 2, . . . , 8. It seems to be the case that the authors of [7,10] 
have not been aware of [8,9]. Some examples of explicit solutions to (1.1) were provided in 
[7,10] under the restriction that the refiection and connection coefficients are all zero. We 
should indicate that the terminology used in [7,10] differs from that used in [8,9]. For the 
benefit of the reader the relationships between the quantities in [7,10] and those in [8,9] 
are indicated in Section 5 of our paper. 

It is already known [5,10,12] that (1.1) is related to the coupled system of nonlinear 
partial differential equations given in (4.2), which is solvable by the inverse scattering 
transform method and is related to the Gelfand-Dickey hierarchy [5]. The time evolution 
of the scattering and bound-state data forms a crucial step in that inverse scattering 
transform. 

One of our goals in this paper is to analyze the bound states for (1.1) by analyzing it in 
the sector of the complex A;-plane, which was also the domain used in [8,9], and to obtain 
the time evolution of the corresponding bound-state data; namely, the time evolution of 
the dependency and norming constants. The time evolution of the bound-state norming 
constants is given in [10] by analyzing (1.1) in the entire complex /c-plane. Another goal 
of ours is to obtain the time evolutions of the scattering and connection coefficients as 
well as other coefficients related to scattering, which is absent in [7,10]. In our paper 
we also clarify how the relevant quantities in [7,10] are related to those of [8,9], and we 
provide some illustrative examples of explicit solutions to (1.1) helping to understand the 
corresponding scattering and bound states better. 

Our paper is organized as follows. In Section 2 we provide the preliminaries and 
introduce the Jost solutions ip^ and ip- as well as the exponential solutions 0+ and 0_ 
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to (1.1), and we present the scattering data in terms of the spatial asymptotics of the 
Jost and exponential solutions. In Section 3 we analyze the bound states associated with 
(1.1) and introduce the dependency constants and norming constants for each bound state. 
In Section 4 we consider the system of integrable nonlinear partial differential equations 
satisfied by the time-evolved potentials u{x, t) and v{x, t), and we obtain the time evolution 
of the corresponding scattering and bound-state data. In Section 5 we explain how the 
quantities associated with (1.1) and used in [7,10] are related to those used in [8,9]. Finally, 
in Section 6 we provide some illustrative explicit solutions to (1.1). 

2. PRELIMINARIES 

Consider the case that the potentials u and v in (1.1) are functions of x alone and 
do not depend on the parameter t. In Section 4 we will consider the case where u and v 
depend on both x and t. 

As in [8] let us introduce the Jost solutions V'+ and ip- to (1.1) in the sector argA; e 
(0,7r/4) with the asymptotics 

ip±{k,x)^e^'^''[l + o{l)], x^±oo, (2.1) 

in such a way that e^^''^ip±{k, x) remains bounded for all x E H. Similarly, let us in- 
troduce the exponential solutions and 0_ in the sector arg/c e (0,7r/4) satisfying the 
asymptotics 

(l)±{k,x) = e'^'"=[l + o{l)], x^±oo. (2.2) 
It follows that these solutions satisfy the respective integral relations given by 

Mk,x) = e-'^- + wJ^J f^^'""^^^ " ^^'""'^^^ + kjv{y)]e''^^^-yU+{k,y), 
V'+(fc,x) = e^^- + — / dy[e^u{y)-k'^u'{y) + k,v{y)]e''^^--y^,p+{k,y) 

~w£ ^^[^3«(y) - kiu'iy) + k^v{y)]e^-^--y^+{Ky), 
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iP-{k,x) = e-''''= + 



Ak^ 



dy [kfu{y) - klu'iy) + k^v{y)]e^'^^-y^_{k,y) 



^ £ dy[e^u{y)-k]u\y) + k^v{y)]e^^^^-y^_{k,y), 



i=2,3,4 



where we have defined 



k\ := k, k2 '■= ik, ks := —k, k^i := —ik. 



(2.3) 



The spatial asymptotics of the Jost solutions and the exponential solutions can be 
obtained with the help of the four integral relations given above. Below, we list those 
asymptotics, through which we introduce the scattering and connection coefficients for 
(1.1). For argA; e (0,7r/4) we have 

( ip±{k,x)^e^'''''[l + o{l)], x^±oo, 

1 



il^±(k, x) = e 



T{k) 



+ 0(1) 



X — > ^oo. 



(/)±(A;,a;) = eT'=^[l + o(l)], x ^ ±oo, 

(l)±{k,x) = e'^''''[A{k) + o{l)], x ^ ^oo. 
For arg A; = we have 

( ip±(k,x) = e'^''''^[l + o(l)], x^ioo. 



i^±{k, x) = e 



„±ifca; 



1 



+ 



R±ik) 



T{k) T{k) 
^>±(/c,x) = eT'=^[l + o(l)], x^±oo, 
{ (j)± (k, x) = e^^"" [A{k) + o(l)] , a; ^ TOO. 
For arg A; = 7r/4 we have 

{ '4^±{k,x)^ e^'^"" [1 + C± [k)e'fkxTikx + ^(1)] , x^ ±oo 

1 



X — > TOO, 



ip±{kj x) = e 



„±ikx 



T{k) 



X — > TOO, 



0±(/c,a;) = eT'=^[l + o(l)], a; ^ ±oo, 
I 0±(A;, a;) = e^^^^ [A{k) + 5±(A;)e±'=^±*'=^ + o(l)] , 
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X — > TOO- 



(2.4) 



(2.5) 



(2.6) 



We emphasize that A{k) and T{k) are defined for argA; e [0,7r/4], R+{k) and R-{k) 
are defined for argfc = 0, and the four coefficients B^{k), C+{k), and C-{k) are 

defined only for arg /c = 7r/4. The coefficient T is known as the transmission coefficient, i?^. 
and R- are the left and right reflection coefficients, respectively, C+ and C_ are known as 
the connection coefficients, and 5+ and S_ are some coefficients that can be expressed in 
terms of A, C_|_, and C_, as we will see. It is either known [8] or can easily be shown that 
for arg /c = we have 

R±{k) _ R^{kY 



1 + 



|i?±(/c)p 



1 



\T{k)? \T{kr 
and for arg /c = 7r/4 we have 



C-{k) 



T{ky 



iT{k)*c+{ky 



A{k) = A{ky 



Tik) 



C-{k)\^\C+{k)\, 
A{k)* + B±{k)C±{k), B±{k)* + A{k)C±{k) = 0, 



\B.{k)\^\B+{k)\ 
1 



^ ^ = \A{k)f - \B^{k)f = \A{k)\' (1 - \C^{k)f) , 



where the asterisk denotes complex conjugation. 



(2.7) 



(2i 



3. BOUND STATES 

Since the coefficient of the third derivative in (1.1) is zero, it follows from the general 
theory of ordinary differential equations that the Wronskian of any four solutions to (1.1) 
is independent of x, and that Wronskian is zero if and only if those four solutions are 
linearly dependent. Recall that a Wronskian is deflned with the help of a determinant. 
For example, the Wronskian involving the Jost and exponential solutions is given by 

^+ ^- <i>+ (t>- 



W4[V'+,V'->+, 



V'V V'- (t>'+ 4>'- 
n i^'L K 0- 
V'T V'- 0+' 0- 
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Using (2.4)-(2.6) we obtain 

W4iP+{k,x),iP_{k,x),4>+{k,x),4>-ik,x)] = -IQik'^P^, arg/c e [0,7r/4]. (3.1) 

The linear independence and boundedness properties of various solutions to (1.1) help 
to identify bound-state solutions. Recall that eigenfunctions of C correspond to square- 
integrable solutions to (1.1), which are also known as bound-state solutions. It is easy 
to verify that C is selfadjoint and hence its eigenvalues can occur only for real values of 
k"^, i.e. when k lies on the boundary of the region Q introduced in Section 1. Thus, any 
positive eigenvalue of C can occur only on the ray argA; = and any negative eigenvalue 
can occur only on the ray arg k = iv / 4. If A{k) = Oat some fc- value on the boundary of the 
region O and a square-integrable solution to (1.1) at that /c- value does not exist, then we 
call that /c- value a spectral singularity of (1.1). U A{k) = at some /c- value in the interior 
of f2, then we call that A;- value a nonspectral singularity of (1.1). By a singularity we refer 
to either a spectral or nonspectral singularity. It is already known that at a singular point 
the two integral relations given in Section 2 for '^+(^7 x) and i/j-ik, x), respectively, are not 
solvable [8]. Spectral and nonspectral singularities for (1.1) may exist, and some explicit 
examples are illustrated in Section 6. 

Note that a bound state in the region can occur only when A{k)/T{k) = some- 
where on the ray arg A; = or arg A; = 7r/4. Otherwise, as seen from (3.1), the four solutions 

V'-, </>+, and (f)- are linearly independent, and the asymptotics of those four solutions 
given in (2.5) and (2.6) indicate that no linear combination of them can decay simultane- 
ously both as a; ^ -|-oo and x — > — oo. 

If there is a bound state at = « on the ray arg A; = 0, then we must have A{k,) = 0. 
This follows from the first identity in (2.7) implying that 1/T{k) ^ and the fact that 
^(k)/T(k) = at the bound state. 

Since four linearly independent solutions to (1.1) must have respective asymptotics 
proportional to e'^*[l -|-o(l)], e~'^^[l-|-o(l)], e''^^[l -|-o(l)], e~''^^[l -|-o(l)] as a; ^ -|-oo, and 
appropriately similar asymptotics as a; ^ — oo, it follows that a bound state at A; = k must 
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decay exponentially as both as a; — > +00 and x — > —00. In this case, we see from (2.5) 
that s'-iid x) are two linearly independent solutions to (1.1) and they do not 

vanish simultaneously both as x ^ +00 and x —00. Thus, from (2.5) we conclude that 
a bound-state eigenfunction </?(«, x) must be in the form 

(f{K, x) = di(f)+{K, x) = d2(f)-{K, x), (3-2) 

for some nonzero constants di and d2. Since any constant multiple of an eigenfunction is 
still an eigenfunction of JC, only the ratio d2/di is relevant and we can call it a dependency 
constant ai k = k, i.e. 

:= 14^. (3.3) 

(p-{K,X) 

where rj{K) is the dependency constant at A; = k on the ray argA; = 0. Defining the bound- 
state norming constants d+ and d- as 



d±{K) 



dx 4>±{k, x)(P±{k, x) 



n-1/2 



(3.4) 



we see that d- = 'r]d+ and that d±{K) 4)±{k, x) is a normalized bound-state eigenfunction 
of the operator C. 

Having clarified the status of bound states on the ray argA; = 0, let us now consider 
the bound states on the ray arg /c = tv/4. If there is a bound state at A; = k on the ray 
arg/c — tt/4:, then we have three possibilities: 

(i) A{k,) = and 1/T{k) 7^ 0. In this case, an argument similar to the case given on the 
ray arg/c = shows that </>+(«;, x) and 4>-{k, x) are linearly dependent and the bound 
state is simple and has the form given in (3.2). Recall that a bound state occurring 
at A; = K is simple if there is only one linearly independent square-integrable solution 
to (1.1) when k = k. 

(ii) A{k,) 7^ and 1/T{k) = 0. In this case, a similar argument indicates that a bound-state 
eigenfunction must have the form 

(P{k,x) = Ci'lp+{K,,x) = C2lp-{K,x), 



for some nonzero constants ci and C2. Since a bound-state eigenfunction is defined up 
to a constant multiple, only the ratio c^jcx is relevant and we can call that ratio a 
dependency constant at /c = k, i.e. 

V^+(«,a;) 

where 7(«:) is the dependency constant a,t k — k on the ray axgk — n/A. Defining the 
bound-state norming constants c+ and c_ as 



c±{k) 



/oo 
dx iP±{k, x)iP±{k,, xY 
-oo 



n-1/2 



(3.6) 



we see that c_ = 70+ and that c±{k) iIj±{k, x) is a normalized bound-state eigenfunc- 
tion of the operator C In this case there is only one linearly independent bound-state 
eigenfunction ai, k = k, and hence the bound state is simple. 

(iii) A{k,) = and 1/T{k) = 0. In this case, 'ip+{K,,x) and '0_(k, x) are linearly dependent 
solutions decaying exponentially as a; — > -|-oo as well as a; — > —00, and (f)+{K,x) and 
4>-{k, x) are also two linearly dependent solutions decaying exponentially as x ^ -|-oo 
as well as x — > —00. On the other hand, 4!+{k,,x) and (f)^{K,x) are linearly indepen- 
dent, and hence they form a basis for the two-dimensional eigenspace consisting of 
bound-state solutions at A; = k. In other words, the multiplicity of the bound state 
in this case is two. The corresponding dependency constants ?](«:) and ■^{k) and the 
norming constants c+{k), C-{k), d+{K), d-^n) are defined as in (3.3)-(3.6). 

4. TIME EVOLUTION OF SCATTERING AND OTHER COEFFICIENTS 

In the previous sections we have considered the Jost and exponential solutions, the 
scattering and connection coefiicients, and the bound states associated with (1.1) in the 
case where the potentials u and v are functions of x alone. Let us now assume that the 
potentials u and v appearing in (1.1) also depend on the extra parameter t. In that case all 
the relevant quantities associated with (1.1) may also depend on t as well. In this section 
we analyze such a dependence on t by interpreting t as the time variable. 
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Now consider the time evolution of the potentials u{x,t) and v{x,t) from their ini- 
tial values u{XjO) and v{x,0), respectively, so that the time-evolved linear operator C 
corresponding to (1.2) is given by 

C:=dt + dMx.t)dx + v{x,t), (4.1) 

where dx = d/dx, and let [5,10,12] 

A := —Sdx — 6u{x,t)dx — 3ux{x,t), 

so that C and A form a Lax pair. As easily verified, the differential operator Ct + CA — AC 
reduces to a scalar multiplication operator, and in fact we get 

Ct+CA-AC = 0, 

which is equivalent to the system of nonlinear evolution equations 

Ut = lOUxxx + QuUx — 24fa;, 

(4.2) 

— ^'^XXXXX ^ ^'^'^XXX ^ ^'^x'^XX QUjU ^ XXX * 

Note that we use subscripts to denote the appropriate partial derivatives. 

Since u and v vanish as x —>■ ±oo at each fixed t, we have A —8d^ as x —>■ ±oo. 
From (4.1) we also see that 

Ct = utdl + Uxtdx + Vf 

In order to solve the initial-value problem related to (4.2), i.e. to determine u{x,t) and 
v{x, t) that solve (4.2) when u{x, 0) and v{x, 0) are specified, we are interested in analyzing 
the time evolutions of the scattering and other coefficients associated with (1.1). 

Towards our goal, we first analyze the time evolutions of the Jost solutions ip+{k, x, t) 
and il)-{k,x,t) and the exponential solutions (f)+{k,x,t) and (f)-{k, x,t), from which the 
time evolutions of other relevant coefficients are easily extracted. 

Theorem 4.1 In the region axgk e [0, 7r/4] the time evolutions of i(^+{k,x,t), il)-{k,x,t), 
(f)+{k,x,t), and (f)-{k, x,t) are given by 

[dt - A]iP± = TSik^^±, [dt-A](f>± = l^8k^(f>±. (4.3) 

10 



Moreover, the time evolutions of various coefficients appearing in (2.4)- (2. 6) are given by 

T{k,t) = T(k,0), A(k,t) = A{k,0), arg/c e [0,7r/4], 

R±{k,t) = R±{k,0)e'^^^'''^\ argA; = 0, 
B±{k,t) = S±(A;,0)e±(s*'='-«^')*, C±{k,t) = C±(A;, 0)eT(«*'='-«'=')*, arg/c = 7r/4. 

PROOF: The proofs in (4.3) can all be given as in the case of the time evolution of iJj+ 
for argA; = 7r/4, which is outlined below. It is known [1-3,11] that [dt — A]ip+ must satisfy 
jC(p = k'^(p, where C is the operator in (4.1). Thus, we have 

[dt - A]ip+ = ci(k,t)ip+ + C2ik,t)ip- + C3ik,t)(f)+ + c^{k,t)(j)-, (4.4) 

for some coefficients Cj{k,t) to be determined. By evaluating (4.4) as a; ^ — oo and 
X — > +00, with the help of (2.6) we obtain 

-jr) -8^/c3^ = Cl^ + C2[e-^'=^ + C_e'=^] + C3[^e-^" + 5+e^^^]+C4e'=^ (4.5) 

(C+)te-'=^-8zA;3e*'=^ - Sk^C+e''''' 

^-ikx (4.6) 

By matching the corresponding coefficients of the exponential terms in (4.5) and (4.6) we 
get 

ci = -8ik^, C2 = C3 = C4 = 0, Tt = 0, {C+)t = {Sk^ - Sik^)C+, 

and hence the first equation in (4.3) for is confirmed when argA; = 7r/4, and we also 
get the time evolutions of T and C+ when argA; = 7r/4, as stated. The remaining parts of 
the proof are obtained in a similar way. | 

The implication of Theorem 4.1 that T{k, t) and A{k, t) do not change in t is signifi- 
cant. As we have seen in Section 3, at a bound state k = k we must have A{k, t)/T[K, t) = 
0, and at a singularity k = Kwe must have A{k, t) = 0. Hence, the /c-values corresponding 
to bound states or singularities of the operator C of (4.1) also remain unchanged in time. 
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Theorem 4.2 Assume that k = k corresponds to a bound state of (1.1). The time evolution 
of the bound-state dependency constants 7(«;, t) and 'q{K, t) and the evolution of the norming 
constants c± {k, t) and d± {n, t) are given by 



c±{n, t) = c±{k, 0) e±4(i+i)«'t, d±{K, t) = d±(«, 0) e±4(i+*)«'*, 
7(«, t) = 7(^v, 0) e-8(i+*)«'*, ?7(«, t) = ri{K, 0) e'^^^+'^^'K 



(4.7) 
(4.8) 



PROOF: Let us assume that there is a bound state at k = k, somewhere on argfc = 7r/4 
with 1/T{k,0) = 0. Then, 'iIj^{k,, x,t) is a bound-state solution and the norming constant 
c+{K,t) can be defined as in (3.6) via 

-1/2 



c+{K,t) := 



dx iP+{k,, X, t)'ip+{K,, X, ty 



(4.9) 



so that c_|_(«:, t)'0_|_(K, x, t) is normahzed, i.e. its L^-norm is equal to one. Let us now find 
the time evolution of c+(k, t). From (4.3) and its complex conjugate we obtain 



[df + Sd^ + 6dxu{x, t) + 3ux{x, t)]tp^{K, x, t) — —SiK^'ilj^{K, x, t), 



[dt + 8d'^ + 3dxu{x, t) + 3ux{x, t)]i;+{K, x, t)" = 8i(K*)3^+(^v, x, t)*, 



(4.10) 



(4.11) 



where we recall that the potentials u and v are assumed to be real valued. Multiplying 
(4.10) by x,t)* and (4.11) by x,t), and adding the resulting equations we 

obtain 

dt\ij+f + dx [8r+{dli;+) + 8i;+{dlr+) - 8{dxi^+){dxr+) + Qu\i;+\^] = -8(1 + i)K'\ij+f, 

where we have used the fact that k* = —ik on the ray arg/c — 7t/4. Integrating over the 
real axis and using the vanishing of '0+, dxij+, d^ij+j and u as x ^ +oo and x —>■ — oo, we 
obtain 

t)|2 = -8(l + z)«3 / dx\iP+{n,x,t)\''. (4.12) 

J — oo J —oo 

Using (4.9) we can write (4.12) as 



d 
~dt 



c+{K,ty 



-8(1 + i)K^ 

C+(K,t)2 
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or equivalently we obtain 



4:{l + i)K^C+{K,t), 



dt 



which yields 



C+{K,t)^C+{K, 0) 6^^^+^) 



The time evolution of the norming constants c_(«;,t), and d-{K,t) appearing in 

the analogs of (3.4) and (3.6) can be obtained in a similar way. With the help of (4.3) 
we obtain (4.7), and hence the dependency constants j{K,t) and r]{K,t) appearing in the 
analogs of (3.3) and (3.5), respectively, evolve according to (4.8). | 



As seen from (1.1), if f{k,x) is a solution to (1.1), so are f{—k,x), f{ik,x), and 
f{—ik, x). Thus, a solution to (1.1) known in the region Q in the complex fc-plane can be 
extended to the three regions obtained by rotating Q by n/2, n, and 37r/2, respectively, 
around the origin of the complex A;-plane. Moreover, since the potentials u and v are real 
valued, f{k*, x)* is also a solution and hence a solution known in a region in the complex 
/c-plane can be extended to the symmetric region with respect to the real axis. Thus, 
solutions known in Q can be extended to the entire complex /c-plane. 

In [10], some four solutions ^j{k,x) to (1.1) for j — 1,2,3,4 are presented on the 
whole complex /c-plane with spatial asymptotics 



totics as a: ^ ±oo, we see that those four solutions are related to the Jost and exponential 
solutions appearing in (2.1) and (2.2) as follows: 



5. A COMPARISON OF REFERENCES [8] AND [10] 




— oo. 



where the aj{k) are certain coefficients and the kj are as in (2.3). Comparing the asymp- 



'^i{k,x) 



(t>-{k,x) 

m 



*2(A;,a;) = 'ilj+{k,x), 
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^3{k,x) = (l>+{k,x), ^4{k,x) = T{k)il;_{k,x), 

where A and T are the coefficients appearing in some of (2.4)-(2.6). Then, as k moves to 
the boundary of argA; e (0,7r/4) from the interior, we see that the reflection coefficients 
''o(^)) ''i(^)) and r2{k) defined in [10] are related as follows to the coefficients used in [8] 
and in our paper: 

ro(fc) = R-{k), r,{k) = C+iky, r,{k) = 

Moreover, the quantities aj{k) appearing in [10] are related to the quantities used in [8] 
and in our paper as 

1 1 
'''^^^ " ' ""^^^^ " T(^' ""^^^^ " ""^^^^ " 

Thus, the inverse problem has been analyzed in [10] in the special case R±{k) = C±{k) = 

B±{k) = 0. In that case, A{k) and T{k) simply become rational functions of k with 

asymptotics 1 + 0{l/k) as /c — oo and with appropriate jump conditions on the rays 

arg A; — {I — 1)77/4 for / = 1, 2, . . . , 8. One can then formulate the inverse problem on the 

entire complex A;-plane as a Riemann-Hilbert problem and solve it explicitly. 

6. EXAMPLES 

In this section we present some explicit examples of solutions and relevant quantities 
associated with (1.1). Such examples should help to understand better the scattering and 
bound-state data for (1.1). It is already known [5,8,10] that if f{k,x) is a solution to the 
Schrodinger equation 

-f"ik,x) + qix)fik,x) = efik,x), 
then f{k,x) is also a solution to (1.1) when 

u{x) = —2 q{x), v{x) — q{xY — q"{x), 

because in that case we have 

D^ + DuD + v=(-D^ + qf. (6.1) 
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Note that (6.1) holds in our first and fourth examples below, but it does not hold for our 
second and third examples. 

Example 6.1 Let u{x) = v(a;) = for a; < and 

4e^ 



u{x) 



, , 2e^(l-e^)2 
v{x) = —jr" -TT—, X > 0. 



(l + e^)2' ' ' (l + e^)^ ' 
By using the continuity of il:'^, •i/'" at a; = and the jump condition 

V'^'(A;,0+) - Vr(^,0-) = -^0+) V;(A;,0), 

we can determine all the quantities relevant to (1.1). In terms of 

2i 



f{k,x) := e 



ikx 



1 - 



(2A; + z)(l + e^) 



we have the Jost and exponential solutions 

f{k,x) + C+{k)f{ik,x), x>0, 



ijj+{k,x)^< I 



T{ky T{k) ' 
f{ik, x), a; > 0, 



-ikx 



+ ci{k)e 



kx 



x<0. 



A{k) + B+{k) + C2{k)e^'= + C3{k) e 



-ikx 



x<0, 



where 
A{k)-- 

C2{k) 

C+{k) = 



(8/c2 + l)(2A;-l) 



16fc3 

-2k + 1 



16A;3 



cs{k) 



B+ik) 

(4 + 4i)/c2 -i 
16k^{2k + 1) ' 



(4 - 4i)k'^ + i 
16A;3(2/c + 1) ' 



ci{k) 



{2i - 2)k + {I + i) 



R+ik) 



2k{8k^ + 1) 

16A;4 + (2 + 4i)A;2 - 1 
(8A;2 -i)(16A;4 + 2zA;2 - 1) 



(4 + 4z)A;^ 



T(fc) = 



2k{2k - l){2k + l){2k + i){8P + 1) 



(8/c2 + l)(2A; + i)(2A;-l)' ' ' {Sk"^ - i){16k'^ 2ik^ - 1) 

The remaining coefficients B_, and C_ can easily be evaluated by using (2.7) and (2.8). 
In this example, there is exactly one simple bound state at A; = k with k := (1 + z)/4, 
where T{k) has a simple pole. A corresponding bound-state eigenfunction is a constant 
multiple of ^)) ^ind it can be chosen as 



-x/4 



{cos(a;/4) + sin(a;/4) + [cos(a;/4) - 3 sin(a:/4)]} , x>0, 



(f{K, x) = < 1 + 

e^/^ [cos(a;/4) - 3 sin(a;/4)] , a; < 0. 
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Even though A and T each have a zero at A; = 1/2 on the ray argA; = 0, their ratio A/T is 
nonzero at A; = 1/2, which corresponds to a spectral singularity and not to a bound state. 

Example 6.2 Let u{x) = and v{x) = —e5{x), with 5{x) denoting the Dirac delta 
distribution and e being a real, nonzero parameter. We want i/j" to be continuous 

at x = 0, and ^p"'{k, 0") = -e + ij"'{k, 0+). We find 

4/C + [i , -ikx _ ^ kx 

I 4A;3 + e ^ Ak^ + e Ak^ + e ' - ' 

r e-'^^, a; > 0, 

Mk,x) = I _ ikx _ _^kx , 

I 4A;3'' ^4A;3'' 4k^ ^ Ak^ + e ' 

The coefficients related to the corresponding scattering problem are given by 

1 _ Ak^ + {l-i)e _ ie 

T{k) Ak^ + e ' ^ 4/e3+(l-z)e' 

Ak^^-e n\ n\ " 

Mk) = —rT:^, B±ik) = - — , C±{k) = 



Ak^ ' ' Ak^' ' 4A;3 + e' 

Note that A{k)/T{k) vanishes on the rays arg/c = and arg/c = n/A only when 

Ak^ + {1 - i)e = 0. 

Since we assume e 7^ 0, we find that there are no such k values if e < 0, and there exists 
exactly one k value lying on the ray arg/c — tv/A when e > 0. Denoting that /c- value by k,, 
we obtain a bound state of multiplicity one at A; = k, where 

1 + i 



K := 



■re. 



2 

Thus, a bound-state eigenfunction is obtained as 

{^iKx _j_ Iq-kx^ X >0, 

^-iKx ^^Kx ^ a;<0. 
Since (ix = 4/ a normalized bound-state eigenfunction is given by 

^x/2 [cos(^a;/2) + sin(^a;/2)] , a; > 0, 
6^e^x/2 [cos(^a;/2) - sin(^a;/2)] , x < 0. 
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<f{K,x) = 



Example 6.3 For any positive constant c, consider 

16c2[l + v/2cosh(2ca;)] 



v{x) = 
We then obtain 



^ ^ [V2 + cosh(2cx)]2 ' 

4c^[V2cosh(6ca;) - 12cosh(4ca;) - 5v^cosh(2ca;) + 4] 
[V2 + cosh(2ca;)]4 



[k - {1 + i)c\[k - {1 - i)c\ 



T{k) 



[k + {I + i)c][k - {I - i)c] 



[k + {I + i)c][k + {I - i)cy ' ' [k - {I + i)c][k + {I - i)c] 

R±{k) = 0, C±{k)=0, B±{k)^0, 



^+{k,x) = e^^'' 
where we have defined 



1 + 



+ 



k + {l + i)c k-{l-i)c 



4)+{k,x) = ilj+{ik,x), 



a{x) 



-V2c- (1 



i]ce 



-2cx 



\p2 + cosh(2ca;) 

This example was presented in [10] in different terminology. Note that A{k)/T{k) has a 
double zero at k = k with k := {l + i)c, which corresponds to a bound state of multiplicity 
two. Two linearly independent cigcnfunctions are given by x) and 4'+{k, x), or they 

can be chosen as real valued, e.g. as 

(v^e-'^^ + 26^=^) cos(ca;) + v^e-^^^ sin(ca;) 



(Pi{k,x) = 



a/2 + cosh(2ca;) 

(-\/2e~'^^ + 2e^^) sin(ca;) — A/2e~'^^ cos(ca;) 



V2 + cosh(2ca;) 



Example 6.4 Consider the potentials 



u{x) 



v{x) 



1x1 + 1)2' {\X\ + IY' 

Using (2.4)-(2.6) and the continuity of the solutions to (1.1) and the continuity of their 
first, second, and third x-derivatives, we get 

k{k^ + k^ + k'^ + 2k + 2) 



T{k) = 



{k + i)[k'^ + (1 + i)k^ + ik^ + (1 - i)k + 3] ' 
17 



C+{k) 



Aik) = 

-2{k + i 



k'^ + k^ + k^ + 2k + 2'' 



B+{k) = 



-2i{k + l){k + i) 



i{k^ + k + 3) 



i(^+{k,x) 



{k + i)[k^ + (1 + i)k^ + ik^ + (1 - i)k + 3] ' 
f{k,x) + C+{k)f{ik,x), x>0 

1 f^-k,-x) + ^^f{k,-x)+C+{k)f{ik,-x), x<0, 



T{k) ' ' T{k) 
and for a; > we have 4>+{k, x) — f{ik, x), while for a; < we have 

4>+{k, x) = A{k) f{-ik, -x) + B+{k) f{-k, -x) + B+{ky f{k, -x) 



e-2 
k^ 



f{ik, -x) 



where we have defined 



f{k,x) := e 



ikx 



1 + 



k{x + l) 

In this example there exists exactly one bound state at k = n with k := 0.778(1 + i) 
corresponding to the simple pole of T{k) on the ray arg/c = tv/4. Note that we have used 
an over line to denote the round-off on the digit. An eigenfunction for that simple bound 
state is a constant multiple of ip+{K,, x). At k = 0.476 + 1.183z in the interior of the sector 
Q, both A and T have simple zeros without A/T vanishing there; thus, that A;- value does 
not correspond to a bound state and it corresponds to a nonspectral singularity. 
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